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Abstract. In [13], K. Roth showed that the expected value of the discrepancy of the cychc 
shifts of the N point van der Corput set is bounded by a constant multiple of ^JTogN, thus 
guaranteeing the existence of a shift with asymptotically minimal discrepancy, [11]. In 
the present paper, we construct a specific example of such a shift. 



1. Introduction 

Let c [0, 1]^ be a finite point set of cardinality N. The extent of equidistribution of 
can be measured by the discrepancy function: 

D^^{x^,X2) := tl(^N n [0,xi) X [0,X2)) - Nxj • X2, 

i.e. the difference between the actual and expected number of points of in the rectangle 
[0, Xi) X [0, X2). The main principle of the theory of irregularities of distribution states that the 
size of this function must increase with N. The fundamental results in the subject are: 

K. Roth's Theorem. ([11], 1954) For any set c [0, 1]^, we have 
(1-1) ||D^J|2>(logN)i/2 

where ">" stands for "greater than a constant multiple of". 

W. Schmidt's Theorem. ([14], 1972) For any set J^ljv c [0, 1]^, we have 

(1.2) \\D^J^>\o^N. 



Both theorems are known to be sharp in the order of magnitude (e.g., [6], [7], [12], 
[2]). One of the most famous examples, yielding sharpness of (1.2), is the van der Corput 
"digit-reversing" set, [6]. For N = 2" points, it can be defined as 

(1.3) 'Vn = {(0.flifl2 • • • ^^nl/ 0.fl„fl„_i . . . aiail) : ai = 0,1} , 
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where the coordinates are given in terms of the binary expansion. Unfortunately, most 
"classical" sets with minimal L°° norm of the discrepancy fail to meet the sharp bounds in 
the norm. In fact, Halton and Zaremba [8] proved that 

(1.4) ||D^„||2 = ^ + 0(n)^(logN)2. 



There are three standard remedies in the theory for this shortcoming. To achieve the 
smallest possible order of the discrepancy, one can alter the sets in the following ways: 

1. Davenport's Reflection Principle. Informally, if P has low L°° discrepancy, then the set 
P = PU{{l-x,y) : (x, y) G P} has low discrepancy. This was demonstrated by Davenport 
[7] in the case of the irrational lattice, and by Chen and Skriganov ([5], see also [10]) in the 
case of the van der Corput set. 

2. Digit Scrambling. This procedure, initially introduced in [3], has been extensively 
studied; a comprehensive discussion can be found in [9]. 

3. Cyclic shifts. This transformation is the subject of this paper. It has been proved by 
Roth, [13], (see also [12], where the translation idea was originally used), that for the cyclic 
shifts of the van der Corput set 

(1.5) = [({X + a) mod 1, y) : (x, y) e , 

the expected value of the discrepancy over a satisfies 

(1.6) f\\D^^J\2da<n"' = (\osNf\ 
Jo 

This implies that there exists a particular cyclic shift of the van der Corput distribution with 
minimal norm of the discrepancy function. However, this was purely an existence proof 
and no deterministic examples of such shifts have been constructed. In the present paper, 
we "de-randomize" this result and provide an explicit value of a, which asymptotically 
minimizes ||D^a||2. We prove the following theorem 

1.7. Theorem. For «o = 1 ~ |r/ where k e'N, in the binary form, is given by 

(1.8) := ( 000 . . . 00 00001111 . . . 00001111 000111 . . . 000111 )^ + 1, 

«o digits 112 digits Jij digits 

with no + wi + n2 = n, ^ = y|/ and no < 568, the cyclically shifted van der Corput set "V^" satisfies 

(1.9) \\D^o\l<n'/^ = {logNf\ 
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Remark. The in the end of (1.8) is just a minor nuisance, which simplifies some 
calculations, and is not important. In fact, one can easily see that a cyclic shift by the 
amount a = 1/N = 2"" changes the discrepancy by at most 1 at each point. 

We would like to point out that most constructions of sets with minimal order of U 
discrepancy (which are important in applications to numerical integration) are probabilis- 
tic; explicit constructions are rare. In fact, the first deterministic examples of such sets 
in dimensions d > 3 have only been obtained quite recently by Chen and Skriganov ([4], 
[15]). 

The outline of the paper is the following: in §2 we deal with the quantities j^^ D^^Xx)dx 
and D'y«{x)dx (which can be viewed as the "zero-order" term of the expansion in any 
reasonable orthonormal basis) and minimize the latter. In §3, we examine the Fourier 
coefficients D^,_(ni,n2) when (ni,n2) ^ (0,0) and show that they do not change too much 
under cyclic shifts. 

We will refer to the two parts of the discrepancy function as "linear" and "counting": 

(1.10) Lm{Xi,X2) =NXi-X2, 

(1.11) C^,(Xi,X2)= 

l[pi,l)x[p2,l)(^l'^2)- 

In proving upper bounds for the discrepancy function, one of course needs to capture a 
large cancelation between these two. 



2. The integral of the discrepancy function 



Recall that in our definition of the van der Corput set, 'V,, = {(O.fli . . . a„l, 0.a„ . . . a2a-[l)}, 
both coordinates have I's in the (n -I- 1)^' binary place. This is just a technical modification, 
which ensures that, for any a = j/2", j G Z, the average value of both coordinates in "V^ is 
one-half: 

1 v-- 1 V ^ 

(^•1) 2n / 1 ~ 2" / ' ~ 2' 

This makes many formulas look 'cleaner' and is not essential to the computations. 



It has been noticed (see [8], [1]), that the quantity D^,,(x)dx is the main reason why 
l|D^J|2 is large. Indeed, if one compares (1.4) and (2.4) below, it is easy to see that 



(2.2) 



D. 



I 



D 



'V,, 



We include the proof of the lemma below for the sake of completeness. 
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2.3. Lemma. For the van der Corput set 'Vn 



(2.4) 



/ 

-^[0,1)' 



D^„(^) dx= -. 



Proof. The linear part of the discrepancy function clearly gives us 



(2.5) r Lm = 2"-2 . 

-'[0,1)2 



Let Xi, . . . , X„ be independent random variables taking values {0, 1} with probability j. 
A straightforward computation yields 



I C^„ (xi, X2) dxi dx2 = (1 - Pi)(l - P2) 



= 2"E 



1 - ^ Xj2-' - 2'"-^ 



n+k-l i-f-n-l 



k=l 



(2.6) 



Combining (2.5) and (2.6) proves the lemma. 



□ 



In what follows we prove that the average of j^^ D^y^dx over a is zero. Besides, we 
construct a specific value of ao, for which 



r 



Dn,«o dx X I. 



2.7. Theorem. Assume that a e [0, 1) is an n-digit binary number. Then 



j 

J [0,1 



(2.8) E„ f D^adx = 0. 

U]2 



Proof. We denote 1 - a = ^ (k = 1, . . . , 2") and start with the following computation: 



r Ct.;; = y (1 - Pl)(l - P2) 
J[0,1V pt^a 
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(l-pi-«)-(l-p2) + Yj (2-pi-«)-(l-p2) 

pe'V„:pi<l-a pe'V„:pi>l-a 



j 



Cn.„dx + (1-«)^(1-P2)- Yj 



P'l]^ pe^n pe'Vn-.piKl-a 

(2.9) = r C^,dx - ^ + P2. 

Next, we examine the behavior of the last sum above. Using the structure of the van der 
Corput set, we can write 

n 

(2.10) Y v2 = Y^-'m+^'2--"-\ 

pe^„:pi<k/2" 1=1 

where k2~"~^ comes from the final I's in the expansion of p2 and 

(2.11) fiik) = #{0 < < - 1 : such that the f' (from the end) binary digit of / is 1}. 
It can be seen that 

(2.12) fi{k) = 2^-1 m if - 1 = 2'm, 2'm + 1, ... , 2'm + 2'"^ - 1, and 



(2.13) fi{k) = + j if - 1 = 2'm + 2'"^ + / - 1, where 1 < ; < 2^~\ 

Thus, if we set fi{k) = 2^-^mi{k) + jiik), where < mi{k) < 2""' and 1 < jiik) < 2'"^ we have 
Efcm/(fc) = \ ■ (2"-' - 1) and 

E,;Kfc) = \ ■ \{2'-' + 1), 
where the extra one-half above comes from the fact that jj{k) = half of the time. Thus 

(2.14) Ekfiik) = 2"-2 - 2'-^ + 2'-3 + 1. 
Plugging this into (2.10), we obtain 

pen^„:pi<<:/2" /=1 

(2.15) =2«--rf 

Finally, equation (2.9), together with (2.15) as well as (2.4), yields 

E r = r D^, - e4 + y Pi 

(2.16) =^_(2«-.l).(2-_^.l) = 0. 
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□ 

To facilitate the construction of an example, we further look at the functions fi{k) = 
2'"^m;(fc) + ii{k) defined above, (2.11)-(2.13). Assume that - 1 is written in the binary 
representation: 

n 

k-l=Y^kj-2'-' = (k„K.^...k2h\. 

7=1 

By construction, mi{k) = Qcnkn-i ■ ■ - ki+i)!, besides, when fc/ = 0, we have //(fc) = 0, and if 
ki = 1, ji{k) = . . . ki)2 + 1. Thus, fi{k) can be written in closed form in terms of digits of 
- 1 as follows 

n 1-1 

(2.17) m = kj2i-' + krYj ^i^''' + ^i- 

j=i+i j=i 

Indeed, if ki = 0, the last two terms will disappear, otherwise, theyTl equal exactly ji{k). 
Plugging this into (2.10), we obtain 

n-l n n n /-I 

(2.18) p^=LL^r +L^'-2-^+YuL^rkr 2'-^-'. 

pe^„:pi<k/2" 1=1 i=l+l 1=1 1=2 j=l 

Obviously, the second term above is bounded by one. Next we shall look at the first term 
in (2.18). At this point we assume that 

n 

(2.19) Y,kj = 1 + 0(11 

7=1 

i.e. approximately half of the binary digits of - 1 are ones and half are zeros. We have 

1=1 j=l+l 7=2 1=1 j=2 

(2-20) =\L^i-^'~'-\L^i = 1^-1 + 0(1). 

7=2 7=2 

As to the last term of (2.18), we have the following lemma: 

2.21. Lemma. For every n G N, there exists k : 1 < k < 2" with Yj^^i = + 0(1), where 
k-1 = (k„k„.i . . . ^2^1^^, so that 

n 1-1 

(2.22) YY,k..kr2'-'-' = ^ + 0{l). 

1=2 i=l 
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Assuming this statement for the moment, putting together (2.18), (2.20), and (2.22) for k 
defined by Lemma 2.21 above, we obtain 



Z/1 n\ n In 

pe^n-Pi<k/2" 

and together with (2.9), (2.6), this yields: 

(2.24) C^adx = (2-2 + + - i) + 0(1) = 2"-' + 0(1). 
Finally, (2.24) and (2.5) give 

(2.25) r D^«(x)rfx = 0(1). 

J [0,1? 

Thus, it remains to prove Lemma 2.21. We shall denote 

n ;-i 

(2.26) S{n, k-l):=YjYj^r^^- 

and will look at some base examples first. Let k' be of the form 

(2.27) k' := (000111 . . . 000111)^, 

where the sequence 000111 is repeated n' times, n = 6n'. We then have the following 
calculation: 



n'-l 



I'-l 



(2.28) S{6n',k') = ^ ^ (2-^''^'^ + 2-^"'^^^ + 2-^"'^^^) ■ ^ [2'''^' + 2'''+^ + 2'''+'] 
i'=i 



l'=0 
n'-l 

2- 



1 "'-1 



/'=! 



6/' . (26'' _ 1) . _1_ . (2-1 + 2-2 + 2-3)(2i +2^ + 23) + i • ^n' 



^ + ^).^ + 0(l) 
72 72/ 6 ^ ^ 



(2.30) 



13 
108 



n + 0(1), 



where the term in (2.28) describes the interactions of digits in different triples and (2.29) 
arises from interactions within the triples. (Notice that the obtained fraction ^ ~ 0.12037... 
is quite close to the desired | = 0.125.) 
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Next we setk" = (00001111... .00001111)2, where the string 00001111 is repeated n" times. 
An absolutely analogous computation yields: 

n"-l 

S{8n",k") = IYj ■ " ■ 2^ ■ ^ ^ ^ + 2' + 2' + 2^) 

;"=i 

I'^/l ll 1\ /I 1 1 



2>^V2 \2 4/ \2 4 8 



(2.31) 



19 
136 



n + 0(1). 



We are now ready to define the number k which satisfies (2.22). Set 

(2.32) k-l:=i 00001111 . . . 00001111 000111 . . . 000111 ) . 

' " ^ ' 

ni digits «i digits 

Then we have, 

(2.33) S(ni + n2, - 1) = S(ni, k!) + S(n2, k!') + J(ni, nj), 

where l(ji\,n2) describes the interaction between the two parts of k. We can estimate: 





( n \ 


' "1 




(2.34) J(ni,n2) = ^- 








< 




J=f!l+1 , 









We now choose Ui and n2 so that ^ = f|, i-e. Wi = H'^/ ^2 = ^'^^ We then obtain 



13 19 
S(n,t -!) = -„, + -„, + 0(1) 



13 -54 19 • 17 
+ 



108 • 71 136 • 71 



)n + 0(l) 



(2.35) 



n 



= 3+0(1), 



which finishes the proof of Lemma 2.21. Thus, if we set ao = 1 - |r, where k is as defined 
in (2.32), then the cyclic shift of the Van der Corput set by «o satisfies 

(2.36) r D^,{x)dx = 0{l). 



Remark. Of course, the above construction only works when n is a multiple of 71 • 2 • 4 = 
568. However, it can be easily adjusted for other values of n just by setting the "remainder" 
digits equal to zero. 



CYCLIC SHIFTS OF THE VAN DER CORFUT SET 

3. The Fourier coefficients of the discrepancy function 



9 



Having eliminated the main problem, we shall now proceed to show that the remaining 
part of behaves well under cyclic shifts. We shall use the exponential Fourier basis 
(rather than the more standard in this theory Haar basis) since it is better adapted to cyclic 
shifts. 



Obviously, for any a, we have 



2"-l 



(3.1) ^ g-27impi _ ^ g-2nimpi _ '^^^-Ini^j , g-ni§t _ ^ 
pe'V" pe'Vn j=0 



0, if m ^ mod 2", 
m, ii m- 2"m', m'- even, 
-m, if m = 2"m', m'- odd. 



Fourier coefficients in the case n\, nz 0. We first note that, for n\, rii 4^ 0, the Fourier 
coeficient of the linear part is: 

— N 

(3.2) LN{ni,n2) = . 

The counting part yields 

(3.3) C^Anu m) = - £ (1 - e-^-"") (l - , 
and, thus, 

(3.4) D^Anx, na) = -r^ Y (e'^^^iPi + e-2"^«2P2 _ g-27ii(nipi+„2P2)\ 

We now consider cases: 

• Both n\ and na = mod 2". Then D<v„{ni,n2) = C^:J^, where C takes values -3 
or 1, depending on whether ni/2" and n2/2" are even or odd. 

• ni ^ mod 2", n2 = mod 2". In this case D^„{ni,n2) = -^^^^ • e-™"2/2". 

• nj^O mod 2", ni = mod 2". In this case D^„{ni,n2) = j^^^ • 

• ni, n2 ^ mod 2". Now we have D^„(ni, n2) = -^^^ Lpe<v^ g-z^Xmpi+napz). 

Changing pi to (pi + a) mod 1 in the above computations, with a - j/2", we notice that 

(3.5) |DT/a(ni,n2)| = |D^„(ni,n2)| when ni,n2 4 0. 

Indeed, in the first three cases the coefficient does not change, while in the last it is 
multiplied by e-27J^«i«. 
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Fourier coefficients in the case n2 = 0,n^ 0. We first note that, in this case 

(3.6) Mni,0) = and C^Mi^O) = E " ^"'"'"''O " P^) • 
Thus, taking into account (2.1), we have 

(3.7) d:;: (m, o) = ^ e-^--^^^^^ ■ (i - . 

And, once again, we obtain that 

(3.8) D^; (ni, 0) = D^, (ni, 0) ■ e-^™"^", i.e. |d^; | = |d^J if ni ^ 0, n2 = 0. 
Fourier coefficients in the case rii = 0,n2 0. As above, we can compute 

(3.9) 5:^,(0,^2) = ^ 2], (1 - Pi) ■ 

In the case n2 = mod 2", we obtain, using (2.1), 

(3.10) 5:^,(0,^2) = D:i:;(0'«2) = ^ • 

The only somewhat non- trivial case is when tii - 0, n2 ^ mod 2". The Fourier 
coefficient in this case is 



D^,(0,n2) = ^X(l-pi)-.-^™-''^ 



2nin2 ^ , 

pe'V- 

-1 

(3.11) =D^„(0,n2) + V e-2™"2P2^ where ;c/2" = 1 - a. 

^ pe^„:pi>J:/2" 

We shall examine the last sum above. Assume n2 = 2^m, where < s < n, m is odd. Let us 
look over the part of the sum, ranging over a dyadic interval of length 2~', 1 < I < n. This 
means that the first / digits of pi (and thus, the last / digits of P2) are fixed, and the last n - 1 
(the first n - Z of P2) are allowed to change freely. 

2"-'-l 

(3.12) ^ g-2"m2P2 = g-2Tr;2=m(,„_,+i2-"+'-i+-+<,„2-"+2-"-i) . ^-2nim2-"+'^' j _ 
pe^„:pie[q2->,{q+l)2~') 7=0 

It is easy to see that the last sum equals zero when l + s <n; otherwise, its absolute value is 
at most 2""'. We now split the interval {pi > k/2"} into at most n dyadic intervals of length 
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2"', 1 < I < n. We obtain 

(3.13) I Yj e-^™"^P^\ < Yj 2""' = 2'+^ - 1. 

pe^n-Pi>k/2" l=n-s 

That is, for n2 = 2^m, by (3.11) and (3.13), we have 



(3.14) D^«(0,n2)-D^„(0,n2) 



< 



2nn2 nm 



4. Proof of Theorem 1.7. 



For a function f & ([0, 1]^) and S c Z^, we shall denote by fs the orthogonal projection 
of / onto the span of the Fourier terms with indices in S, i.e. 



(4.1) /s(xi,X2) =^ Yj /("i'"2)e 



2ni(niXi+n2X2) 



{ni,n2)eS 



Due to (3.5), (3.8), and Parseval's identity, we have 



(4-2) V" ■ /Z2\|ni=0} 

Inequality (3.14) yields 



^\l«i=0| 



(4.3) 



n-1 



s=0 m odd 

Thus, we see that ||(D^j2;2\(o,o)||2 indeed does not change much under cyclic shift. The 
inequalities above and (2.2) yield: 

(4-4) \Py:'UJ\. S II(d-v.L«„,„JIz + ('°SNr S (logN)"^ . 
Together with the fact that J Dy^o < 1, (2.36), this finishes the proof: 
(4.5) \\Dy.o\\<(\ogNf\ 
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